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AN EMPIRICAL METHOD FOR DETERMINING 

STATIC DISTRIBUTED AERODYNAMIC LOADS ON AXISYMMETRIC 

MULTISTAGE LAUNCH VEHICLES 

By Ralph J. Muraca 
Langley Research Center 


SUMMARY 


This paper presents an empirical method for determining the distributed 
aerodynamic loads of arbitrary configurations in the linear angle -of -at tack 
range. The method is based upon a compilation of wind-tunnel force and pres- 
sure data. Existing test results for geometric components such as cones, tan- 
gent ogives, spherical segments, cylinders, frustums, and boattails have been 
correlated and reduced to the form of generalized loading functions. These 
generalized loading functions cover a wide range of Mach numbers and geometric 
parameters and can be quickly applied to arbitrary configurations. Distributed 
aerodynamic coefficients for several configurations have been calculated by 
using the data presented herein, and these results are compared with the results 
obtained from wind-tunnel pressure tests of the specific configurations. 


INTRODUCTION 


% 

The effects of body flexibility have been increasingly recognized as impor- 
tant factors in the dynamic loads and stability studies performed on spacecraft 
launch vehicles. The inclusion of body flexibility has introduced the need for 
an adequate description of the distributed aerodynamic loads acting on the 
deformed body. 

The efforts required to obtain distributed data are considerably greater 
than those for obtaining total aerodynamic coefficients . However, the distrib- 
uted aerodynamic data are essential to the adequate analysis of aeroelastic 
divergence, aeroelastic feedback coupling of closed-loop autopilot systems, and 
dynamic loads of flexible structures experiencing atmospheric disturbances. The 
severity of such instabilities and the load magnifications due to flexibility 
justify the extensive efforts required in generating appropriate distributed 
aerodynamics. Unfortunately, the ability to determine distributed aerodynamic 
loadings has not increased at the same rate as the need for such data. 

Three methods are currently available to the aeroelastic ian for the acqui- 
sition of distributed aerodynamic coefficients. These methods can be classified 
as analytical, empirical, and experimental. 



The analytical methods are based upon the assumptions that the partial dif- 
ferential equations governing the flow around a body and the associated boundary 
conditions can be simplified to an extent whereby they would be amenable to 
solution. Detailed derivations of the various analytical methods along with 
their limitations and ranges of applicability are given in references 1 to 7 . 

A second source of aerodynamic data is empirical methods. These methods 
can be derived from both analytical and experimental considerations. Detailed 
explanations of these methods and their limitations are given in references 8 

to 13 . 

The most accurate sources of distributed loads are wind-tunnel tests of the 
specific configuration under consideration. Thece tests constitute the third 
method available to the engineer for determining distributed aerodynamic loads. 
Either pressure -test or force-test models may be used; however, distributed 
data can only be obtained by using pressure models unless the effect of body 
length on the total coefficients can be determined from the force-test results. 
For a detailed description of the distributed normal-force-coefficient slope on 
a body, pressure tests are the preferred method. However, the time and expense 
required to run pressure tests are considerable and this method cannot usually 
be justified in the preliminary stages of vehicle design. 

This report presents an empirical method of determining the distributed 
normal-force-coefficient slope on arbitrary bodies in the linear angle -of -attack 
range. The method is based upon wind-tunnel results presented in references l4 
to 46. These results have been utilized to obtain generalized loading functions 
which have been categorized by geometric components. These components and the 
ranges of geometric parameters for which the loading functions are available are 
shown in table I. 

The empirical method presented herein is not submitted as a replacement 
for any of the other previously mentioned methods which would be applicable for 
a given design or for wind-tunnel testing but is intended to supplement such 
methods. The correlations presented in this paper provide a rapid means of 
determining distributed normal-force-coefficient slopes for a wide range of con- 
figurations and Mach numbers. 

Included as appendixes A and B to this report are the equations which were 
used to reduce the various wind-tunnel results and a numerical example showing 
the application of this method to a typical launch-vehicle configuration com- 
posed of a blunted cone, a cylinder, a cone frustum, and a second cylinder. 
Further comparisons are also made over a range of Mach numbers for a few spe- 
cific configurations for which wind-tunnel data were available. 
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normal-force-coefficient slope, — — , 1/radian 
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local diameter, meters 
reference diameter, meters 

diameter of the nose cone at the tangency point to the nose 
sphere, meters 

diameter of cylinder following cone, meters 
diameter of cylinder following frustum, meters 
normal force, newtons 

fineness ratio of component, — , calibers 

D 0 

length of component, meters 
free -stream Mach number 
static pressure, newtons /meter^ 

dynamic pressure, i pV^, newtons /meter^ 
radius of sphere or radius of ogive, meters 

local radius of cross section at x normal to vehicle axis of 
revolution, d/2 

reference area, meters^ „ 

velocity, meters/sec 

body coordinate of component, x = \i - p Q , meters 

product of the reference area and the distributed normal-force 
coefficient slope, meters/radian 

generalized loading function, l/radian 
angle of attack, radians 



boattail angle, deg 


frustum angle, deg 


0 n nose semivertex angle, deg 

p free-stream density, kilogram/meter? 


M- 




angle of attack compatible with Cp, radians 


generalized loading function, 



1 /radian 


coordinate along the axis of revolution of the assembled 
components, meters 


reference value of \i at the origin of the component under 
consideration, such that, x = \± - \± Q , meters 


Subscripts : 

a afterbody 

f frustum 

n nose 

1 component preceding the component under consideration 

2 component following the component under consideration 

oo free stream 

T tangency point of sphere and cone 

s sphere 


METHOD OF DATA REDUCTION 


The generalized data which have been obtained are graphically displayed in 

d %r > 
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dx 
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this report as the variation of the generalized loading function 
with the component length. In this expression for the loading function, 

represents the product of the reference area and the distributed 
normal-force-coefficient slope and D is the local body diameter. To convert 
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the loading function to distributed normal force all that is required is multi- 
plication by the local diameter, the local angle of attack, and the dynamic 
pressure of the stream. Geometric scale effects are accounted for and the need 
for defining a specific reference area is eliminated by including the reference 
area term as an integral part of the coefficient. This form for the generalized 


loading function 



is most convenient to the aeroelastician since 


it can be readily incorporated into the various analyses mentioned previously. 


The reference material used in compiling the loading functions presented in 
this report were presented in many diverse forms which necessitated specific 
processing to yield the desired results. The data-reduction operations neces- 
sary for reducing force-test and pressure-test results are recorded in appen- 
dix A. Test results obtained for angles of attack up to 5° were included; 
therefore, the generalized loading functions are applicable in this range. It 
should be noted that a wide range of Reynolds numbers is covered in the refer- 
enced data; however, very little information was available from which Reynolds 
number effects on distributed loads could be obtained. In general, flight 
Reynolds numbers will be an order of magnitude higher than wind-tunnel test 
conditions. At present, neither analytical nor experimental methods are avail- 
able for accurate determination of model to full-scale Reynolds number effects 
on distributed loads. It should be noted that a variety of systems of units 
were used in presenting the referenced data. The units used throughout this 
report are the meter, newton, and radian. Geometric angles are given in 
degrees . 


PRESENTATION AND DISCUSSION OF DATA 


Space vehicles may be comprised of various geometrical components such as 
spheres, cones, cylinders, cone frustums, and other bodies of revolution gener- 
ated by given arcs. Based upon these considerations, the loading functions 
applicable to components such as those mentioned previously have been deter- 
mined and are presented for use with multicomponent vehicle configurations. The 
use of generalized data for individual components to determine distributed 
aerodynamic curves for complete configurations is not theoretically proper in 
the subsonic and transonic ranges as a result of the propagation of downstream 
disturbances upstream. Downstream component parameters will therefore affect 
the loading function on upstream components. However, the extent of this effect 
will depend upon the configuration and will be considered subsequently. The 
loading functions are presented in terms of geometrical parameters associated 
with the component under consideration and the other adjacent components. Fig- 
ure 1 shows the relationship between the component geometry and the various 
parameters used in presenting the data. 


Spherical Nose Segments 

Loading functions for spherical noses are presented in figure 2. These 
data are based upon results of references 44 and 15* An analytical expression 
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involving two Mach dependent coefficients was fitted to the pressure distribu- 
tions of reference 44. The resulting empirical expression was used for com- 
puting generalized loading functions by the process of equation (All) in appen- 
dix A. Generalized loading functions obtained from the analytically described 
pressure distributions were then matched with the maxima of comparative data 
from reference 15 in order to ascertain the previously undefined Mach dependent 
coefficients. 

A comparison of the total normal-force coefficient for a spherical segment 
calculated by using Newtonion impact theory with that obtained by integration 
of the loading function of figure 2(b) at Mach 4. 65 indicates that the results 
obtained by using the loading functions are about 15 percent lower than those 
obtained by using impact theory. 


Cones 

The loading functions for cones are shown in figure 5 for various semiver- 
tex angles 0 n and Mach numbers. These data were obtained from reference 16, 
which presents analytical results, and from reference 17, which presents exper- 
imental results. The method of reference 1 6 does not provide data for those 
Mach numbers where the shock wave is detached. However, reference 17 presents 
wind-tunnel data on cones over the Mach number range from 0.50 to 4.06. It was 
found that, in the region where the analytical solution was applicable, the two 
sets of data agreed very well. The effect of blunting is also presented in 
reference 17. In the transonic range, blunting has little effect on 0 $ for 

all values of 0^. As Mach number increases. Cm decreases as a result of 

blunting. The amount of reduction in Cjj is also a function of 0 n . For 

small 0 n the reduction in C^ a is considerable, but as 0 n increases Cu 

for a blunted cone approaches Cm for a sharp cone. 

^ a 

The loading functions of figures 2 and 5 were used to obtain total normal- 
force coefficients for blunted cones at various Mach numbers and nose angles. 
These calculated values of C^ were then compared with the wind-tunnel data 

of reference 17 and good agreement was observed. In general, this method will 
provide a good description of the distributed normal force on a blunted cone; 
however, the loading function on the cone portion of the sphere-cone combina- 
tion is still considered to be constant with cone length as shown in the deriva- 
tion in appendix A. 


Tangent Ogives 

A tangent ogive is a body of revolution which is generated by revolving 
the area enclosed by the intersection of two equal circles having radius R 
about the common chord of the two circles. This configuration may be seen in 
the following sketch which also gives the relationship between f n , R, and D 0 . 
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A tangent-ogive nose shape consists of the portion of the generated body 
of revolution from the sharp point to the maximum diameter D Q and is of 
length l . 

Very few wind-tunnel test results are available from which to determine 
loading functions for tangent ogives . The curves presented in figure 4 were 
obtained from references l8 to 20. These loading functions are based upon pres- 
sure test results. 

Although the available information is limited, a method utilizing the cone 
data of figure 3 can be used to approximate ogive loading functions. This 
method is based upon the assumption that the loading function at a given ogive 
station is well approximated by that of a cone tangent to the ogive at the sta- 
tion. In this manner, a loading function for any ogive can be determined. 

This method tends to predict higher loading on the forward part of the ogive 
and lower loading on the aft portion of the ogive. It should be noted that the 
data obtained in this manner will have a value of zero at the cylinder tangency 
point. To overcome this inadequacy, the value at the tangency point should be 
obtained from the cylinder-following-ogives data of figure 5- The ogive curve 
in the vicinity of the cylinder should be faired to make the two curves tangent 
at the common intercept. 

A comparison of the loading functions for a tangent ogive obtained by using 
the tangent-cone method and those obtained from the ogive data presented in fig- 
ure 4 is shown in the following plot. This method is intended for use where 
loading functions for a specific ogive cannot be obtained from the data of 
figure 4. 
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Cylinder Following Cones and Tangent Ogives 

Figure 5 presents curves applicable to cylinders following cones and tan- 
gent ogives. Both pressure-test and force-test results were used in generating 
the curves of figure 5- In the transonic range the only data presented are 
from reference 21. This reference presents the results of a series of tests in 
which various geometric parameters were varied and the effects on pressure dis- 
tributions were noted. Reference 21 represents the first means by which the 
engineer might predict the variation of local pressure with component parameter 
in the transonic range. 

Some of the references which present force-test data using stacked models 
fail to provide a data point at x/D 0 = 0 (i.e., forebody alone). Whenever 

this omission occurs for cone-cylinder configurations, the value of Cj^ for 

the appropriate cone was obtained from reference 1 6 or 17 • 

In general for M > 1.5 the value of the loading function at x/D 0 = 0 

increases as cone semi vertex angle increases or Mach number decreases. It is 
apparent from figure 5 that the loading functions show less variation with body 
length as Mach number increases and approach an almost constant value at hyper- 
sonic velocities . In terms of total normal force, the data of figure 5 show 
that the contribution of a cylindrical afterbody increases with increasing 0 n 
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at all Mach numbers above Mach 1.0. As Mach number increases the cylinder 
normal force increases for a given 0 n until it reaches a maximum between 
M = 2 and 3* and then decreases gradually as M increases above M = 3- The 
center of pressure of cylinder afterbodies shows a shift rearward as Mach num- 
ber increases. For example, the center-of -pres sure location for a 6-caliber 
afterbody is between 1.2 < x/D 0 < 1.6 at M = 1.5 and- moves aft to 
2.0 < x/d 0 <2.4 at M = 6.86. 

It has been assumed that the effects of downstream components on the 
loading of the bodies preceding them are negligible; consequently, frustum 
parameters for the case of a frustum following a cylinder do not appear on the 
cylinder loading curves . It should be understood, however, that, even in the 
supersonic range, significant subsonic regions can exist in the boundary such 
that disturbances will be propagated upstream. This is particularly true in 
the case of cylinders preceding frustums . Large effects on cylinder loading 
functions due to Mach number and frustum-angle variation have been shown to 
exist by other investigations. These investigations show that the frustum angle 
can radically affect the cylinder loadings through the mechanism of flow sepa- 
ration. The point at which the flow separation takes place is a function of Qf 
and M. The resultant cylinder loading can be drastically Increased in the 
region of flow separation. When these phenomena occur, the loading functions 
on the frustum will also be affected. A list of references dealing with flow 
separation on cone-cylinder frustums and further discussion of this behavior 
are given in the next section. Insufficient data showing these interactions 
precluded an adequate treatment of these phenomena in terms of both cylinder and 
frustum parameters . 


Frustums Following Cone -Cylinders 


Figure 6 presents the loading functions for frustums following cone- 
cylinders. The distributed normal force on frustums following cone-cylinders 
is a function of four parameters: Mach number (M), frustum angle (0f), cylin- 
der fineness ratio and nose semivertex angle (0 n ). As presented in 


figure 6 each set of curves gives the variation of 




with x/D 0 


for various arrangements of the four parameters. For configurations in which 
the cylinder length was greater than 4 calibers, 0 n no longer significantly 
affects the frustum data and therefore is not included as a parameter. It is 
difficult to generalize with regard to the behavior of the loading function for 
the frustum. Besides the number of variables involved, the possibility of flow 
separation occurring on the cylinder cannot be ignored. Flow separation at 
supersonic speeds is related to the boundary -layer shock-wave interaction in the 
vicinity of the cylinder -frustum junction. Reference 22 indicates that flow 
separation increased with increasing Mach number and decreasing Reynolds number. 


Reference 23 also indicates that, as frustum length or frustum angle 
increases, Cjj increases rapidly with Mach number in the transonic range. As 


9 



mentioned previously, this effect is due primarily to flow separation and would 
he reflected in the distributed data as a sudden increase in loading function 
on the frustum. Unfortunately, prediction of these highly nonlinear effects is 
difficult and no attempt has been made to include them in the loading functions 
presented herein. In general, whenever a configuration is encountered where 
the nose and frustum are within 2 calibers of each other and where the frustum 
angle is greater than about 10° it is quite possible that flow separation will 
be encountered. Reference 19 also shows that, as the ratio of frustum base 
area to cylinder base area increases, the point of flow separation moves forward 
on the cylinder. Other factors also affect separation; see, for example, ref- 
erences 47 and 48. 

The parameters having the greatest effect on the load functions on a frus- 
tum are 0^ and M. The next most important parameter is fa, 1* If f a i is 

great enough so that pressure following the expansion around the cone-cylinder 
shoulder has recovered, a load buildup will occur on the forward portion of the 
frustum. 

The general trend of the frustum loading function is to approach the cone 
value asymptotically as frustum length increases. The loading-function varia- 
tions from those of the cone value usually occur on the forward portion of the 
frustum within 2 calibers of the forebody-frustum junction. 

Included in figure 7 are cross plots of the data of figure 6. These curves 
show the variation of frustum loading functions as f a is varied and all 

other parameters are held fixed. The curve for f a ^ = 0 was obtained from 
cone data. In addition to references 22 and 23, the data of figure 6 were 

obtained from references 12, 21, 24 to 29, and 40. Additional information con- 
cerning the effects of cone angle, cylinder length, and frustum angle on flow 
separation and the resultant nonlinear behavior of Cm can also be obtained 

from these references. 


Cylinders Following Frustums 

The loading functions for cylinders following frustums are presented in 
figure 8. The available data on cylinders following frustums are very limited. 
References 21, 26, 27, and 30 were used to obtain these curves. The data of 
references 21, 26, and 27 were obtained from pressure models. The data of ref- 
erence 30 were obtained from force-test models with stacked sections. 

The parameters which affect the loading on a cylinder following a frustum 
are M, 0f, f a , 1* 9 n , ^l/^o^ and x/d 0 . However, it has been pointed out 

in reference 28 that, if fa,l is greater than about k.5 calibers, the nose 
semivertex angle 0 n no longer significantly affects the loading on the cylin- 
der following the frustum. Consequently, for data obtained from models with 
f a ^q > 0 n and f a ^i do not appear as parameters. The general trends of 

the loading function of a cylinder following a frustum are similar to those of 
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a cylinder following a cone. The parameters having the greatest effect on the 
functions are frustum angle 0^ and the cylinder -diameter ratio ^l/^o* ^ n 

the supersonic Mach number range as Dj_^D 0 increases, the influence of 0-f* 
decreases and, at D-j_/D 0 > 0.75, 9f has no effect on the cylinder loading 

functions. In the event that data for a configuration with a diameter ratio 
less than 0.25 are required, the data for cylinders following cones can be used, 
giving a curve for D]_/D 0 = 0. These data presented in figure 9 represent the 

loading function for cylinders following frustums and are cross plots of the 
data of figure 8, with the diameter ratios Dj_/D 0 varied for specific values 
of 0f and Mach number. 


Boat tails and Rearward-Facing Step 

In the determination of the loading functions for boattail and rearward- 
facing-step configurations, the problems of flow separation, Reynolds number 
effects, and component interactions once again present themselves. The data 
presented in figures 10 and 11 can be considered representative of the load 
distribution on these components. Although interpolation of these data is per- 
missible, extrapolations beyond the limits of the parameters given on the indi- 
vidual curves should be avoided. Potentially large errors can be introduced 
when these data are extrapolated. 

The boattail loading functions are presented in figures 10 and 11. The 
boattail data can be broken up into two parts: boattails following short cone- 
cylinders and boattails following long cylinders. Figure 10 presents data of 
the first type; the parameters of interest are 0-^, M, 0 n , ^a,l- 

In general the boattail loading functions tend to approach positive values as 
boattail length Increases. This trend might be due to the presence of the com- 
ponent following the boattail. This same trend is also seen in the curves of 
figure 11 which presents data for boattails following long cylinders with the 
end of the boattail representing the base of the models. These models were 
sting mounted. Consequently, the sting might have had the same effect on the 
boattail loads as a cylinder. In figure 11 the only body parameter is 0-^ 

since the cylinders used in these tests were long enough that 0 n had no 
effect on the boattail loading. The data of figures 10 and 11 were obtained 
from references 21 and 32 to 3^. 

Figure 12 presents the effect of a rearward-facing step on the loading 
function of a cylinder following a cone. These data were obtained with a model 
in which the step was located 2.4- calibers downstream from the cone -cylinder 
shoulder. The parameters which affect the loading function are Di/D 0 , M, 
f a ^i, and 0 n . However, if the shoulder Is located far enough downstream of 
the cone, the effect of a variation in 0 n will be secondary. Reference b6 

presents the only data which could be used to show the effect of a rearward- 
facing step on cylinder loading functions, consequently, it was not possible to 
determine the effect of a variation of D^/Dq, 0 n , or f a on the loading 

function. 
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Comparison of Wind-Tunnel and Empirical Results 

In an effort to provide some indication of the accuracy of the empirical 
method outlined herein, some comparisons with data obtained from wind-tunnel 
tests of specific configurations have been made. Configurations for which 
wind-tunnel results were classified, or which were used to generate the empir- 
ical loading functions of this report were avoided. 

Figure 13 shows a comparison of the distributed normal-force-coefficient 

f 

slope S S r obtained from wind-tunnel pressure-distribution tests of the 

l ** J 

specific models indicated with similar data obtained through the use of the 
method outlined herein. The configuration shown in figure 13(a) has an upper- 
stage cylinder length of over 4 calibers, and a frustum angle of 5°* For this 
model, flow separation does not occur to a degree where it would drastically 
affect the load distributions. The calculated distributions are in good agree- 
ment with the experimental data except at a Mach number of 1.6. No explanation 
can be given at this time for the poor agreement at Mach 1.6. Some differences 
are also seen for the cylinders just behind the cone and frustum at Mach 2.0 
and 2 . 36 . 

The models In figures 13 (b) end 13(c) have the frustum and nose in close 
proximity to one another, and their respective angles are such that flow sepa- 
ration does occur on the intermediate cylinders just forward of the frustum and 
on the forward portion of the frustum. Along these segments of the model, the 
estimated data tend to predict considerably lower load distributions. For 
those configurations where significant regions of separated flow are likely to 
occur, the predicted loads in general will tend to show fluctuations with less 
extreme peak values. It should be noted that the geometries and Mach numbers 
of figures 13 (b) and 13 (c) represent a severe test for the empirical data of 
this report. 


CONCLUDING REMARKS 


The empirical method for determining distributed loads presented herein 
will yield acceptable results for most configurations and Mach numbers . Care 
must be exercised in applying this method to configurations where flow separa- 
tion is anticipated. The assumption of linearity limits the use of this method 
to angles of attack less than about 5°* The equations used to reduce the basic 
experimental data are included and will allow the user to incorporate addi- 
tional data which might become available. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., September 20, 


1965. 
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APPENDIX A 


DETERMINATION OF DISTRIBUTED AERODYNAMIC COEFFICIENTS 
FROM FORCE- AND PRESSURE-TEST DATA 


Force-Test Data 

Loading functions can be derived from force-test data in which the varia- 
tion of the normal-force coefficient with model length has been determined. 

The following sketch illustrates a typical curve of Cm plotted against x/d o 

CG 

from such a test . 



<3 I 1 


The change in model length is usually achieved by changing the length of 
one of the geometric components used in the configuration through the use of 
stacked segments. The change in Cjj due to the addition of a segment pro- 

CG 

vides a means of determining the contribution of the segment. 

Unfortunately, the curves of plotted against component length were 

not always smooth and determining the slopes of the curves by differentiation 
allowed for potentially large errors. Both graphical and numerical means were 
employed to determine the slopes and the results were weighted in light of con- 
tinuity of data. Where possible, the method of least squares was used to fit a 
curve to the data points and the derivative of this curve was calculated and 
used. The slopes were also calculated by differentiation of the Lagrangian 
interpolation formula. The coefficients developed in reference l4 were used in 
this procedure. The slopes calculated by these various techniques were com- 
pared wherever questionable results were obtained. For those cases where the 
data scatter was especially bad, intuitive estimates were made based upon simi- 
lar results. The exact procedures used for the various component parts of a 
vehicle configuration are subsequently developed. 
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APPENDIX A 


The total normal force is usually given ty 


f N “ 


(Al) 


Taking the partial derivative of equation (Al) with respect to x (where a 
is assumed to he constant) along the body yields. 


dF, 


N 


(dCw 


Sx =qa U 


a 


S + C- 


K 


bS 


a Sx 


(A2) 


The reference area S is allowed to be a possible function of x at this 
point . 


by: 


By definition, the dimensionless loading function, A, is related to F, 


N 


dF; 

dx 


N 


= qa 


{■£} 


= qaAD 


(A5) 


Substituting equation (A3) into equation (A2) yields 

"'ski +s 5r) 


(A4) 


Using the form in which the referenced force-test data for the various com- 
ponents was usually presented, the following specific data-reduction equations 
have been derived. 


Cones .- The total normal-force coefficient Cjj a for cones is usually pre- 
sented as a function of cone semivertex angle and Mach number only. By allowing 
the reference area to be the cone base area, becomes a constant with cone 

ac Na 

length. In this case — ^ — = 0 and equation (a 4) becomes 

ox 


but 


A = 


1 n ds 

5 N aS 


(A5) 


5s _ ]t x 
Sx 2 f n : 


(A6) 


Substituting equation (A6) into equation (A5) yields 
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APPENDIX A 


Cm , 

A = — — - where f = £ = — 
2 f n n D D 0 


(AT) 


Cylinders . - For cylindrical elements, the reference area S is taken as a 
constant equal to the cylinder cross-sectional area. Equation (A4) then becomes 


A £ 

4 a n 


(A8) 


Frustums and boattails .- For both frustums and boattails, the reference 
area is usually taken to be a constant value; therefore, equation (A^) becomes, 
for frustums. 




SCm 


a 


1 + 2 ^ tan 0 f la 


D, 


(A9) 


where D = D Q ■+- 2x tan 0^, and for boattails, 

1 


A = — / 


8Cm 


a 


tan 0-u — 

^o 7 D o 


(A10) 


where D = D Q - 2x tan 0-^. 


Equations (A7) to (A10) were used extensively in reducing the referenced 
force-test data to the loading functions presented in this paper. 


Pressure-Test Data 


The conversion from pressure-test data to the previously defined loading 
function is presented in reference 15 . The final form of the equation relating 
load functions to pressure coefficients as presented in reference 15 is 


1 

D 





0 d0 


(All) 


where 


C P 

0 


test angle of attack 

p - p w 

pressure coefficient, — — °° 

angle between the lift plane and an inclined plane intersecting the 
lift plane along the vehicle axis of revolution 


D local diameter of the cross section normal to the axis of revolution, 

D = 2r. 
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Spherical Segments 

The loading function can be approximated for a spherical segment by use of 
equation (All) and an approximate relationship for the pressure distribution. 

In the following sketchy consider the point b on the spherical surface 
lying in the plane Oab which is inclined the angle 0 with respect to the 
vertical plane and intercepts the vertical plane along the vehicle axis of 
revolution. 




Defining 0 as the angle between the normal to the surface at the stag- 
nation point (ca) and the radius (ab = R) from the origin of the spherical sur- 
face to the point b, then the pressure coefficient Cp can be approximated by 
the relationship 


C p = A cos 2 0 + B 


(A12) 


If a transformation to the x, r, and 0 body coordinates is made, equa- 
tion (A12) becomes 


C-n — A 


l 2 

^1 - — jcos - ~ sin ouf. cos 01 + B 


(A13) 


where 

R sphere radius; R = D 0 /2 

and 

r local radius of the cross section at x normal to the vehicle axis of 

revolution; r = d/2 

x coordinate along the axis of revolution of the vehicle and having its 

origin at the surface of the sphere 
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If equation (A13) is substituted into equation (All) there will result 


cos 0 d0 (Al4) 


I 

II 

§ 

r A [( 

I - l^cos 04 . 

- £ sin cos ( 

(j + B 


' o 


“t 



For the spherical surface; r, x^, and R are constrained hy 




(A 15 ) 


Imposing the constraint of equation (A15) on equation (Al4), integrating, and 
using the customary approximations for small angles of attack (i.e., cos a = 1 , 
sin a = a), equation (Al4) yields 


= ^(1 - tVtRT) 


(ai 6 ) 


where the constant B vanished by virtue of 

/ B cos 0 d0 = 0 
J 0 

Values for the undefined coefficient A are obtainable by comparing the maxima 
of equation (Al 6 ) with comparative maxima obtained from reference 15 for vari- 
ous Mach numbers . 

Equation (Al 6 ) is valid within the bounds on positive x values given by 
the following condition 


x < 1 " s in 9 n 
Do = 2 


(A17) 


where 0 n is the semivertex angle of the nose cone. The preceding limit 

restricts the use of equation (Al 6 ) to x values from zero to the x position 
of the tangent of the nose cone to the spherical surface. 
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APPLICATION OF METHOD 


The empirical method presented in the Body of the report is utilized to 
determine the distributed aerodynamic loads for a typical launch-vehicle con- 
figuration. The configuration for which data are generated is shown in fig- 
ure Bl. Data will be developed only for Mach 2.4 whereas in actual practice 
distributed loads for a range of Mach numbers would be required. 

The vehicle is first broken down into components for which data are avail- 
able. For the configuration under consideration, there will be five components 

(1) Spherical nose cap 

(2) Cone following a spherical nose cap 

(3) Cylinder following a cone 

(4) Frustum following a cone-cylinder combination 

( 5 ) Cylinder following a frustum 


Spherical Nose Cap 


For a spherical nose cap, the bluntness ratio of the conical nose is given 
as b n = 0.40; the cone semi vertex angle is given as 0 n = 22.5°; and the cylin- 
der diameter Dj_ is 2.17 meters. By using the relationships given in table Bl, 


the local diameter at the various x/d q stations can be determined. It should 
be noted that x/D 0 = 0 corresponds to the tip of the vehicle and x/d q = O.3O9 
corresponds to the point at which the sphere meets the cone. For this com- 

The appropriate curve from figure 2(b) is 

values at the x/d q stations of inter- 


ponent, D 0 is the sphere diameter. 


then used to determine the 


{■$} 


est. In this example, interpolation on Mach number is required. These values 
are then multiplied by the local diameters to yield the product of the reference 

area and the distributed normal-force-coefficient slope 
in -table Bl. 


dC w 


'a 


dx 


tabulated 


Cone Following Spherical Nose Cap 

The cone data in figure 3(a) are used to determine the load distribution on 
a blunted cone. Since the loading on a cone is a constant, only two values of 
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are required: one at the tangency point and one at the cone- 

cylinder junction. A straight line is then drawn "between them. The expression 
given in table B 2 is used to calculate the diameter at the sphere -cone tangency 
point. For a sharp cone the diameter would be zero. The reference diameter 
for cones D n is taken to be the maximum diameter of the cone or 

dC %' 
dx 

figure 3 (a)- For this case interpolation is not required on cone semi vertex 
angle. This value is then multiplied by the diameters at the tangency point, 
Drp, and the cone base, D Q , respectively, to obtain the two end points of the 


D 0 = Dp = 2.17 meters. The appropriate value of 


r 01 

4 


is obtained from 


final curve of 


{■¥) 


for the cone. The final curve for the complete 


blunted cone is shown in the following sketch: 




Sphere-Cone Combination 

The curves of figure 5(<i) are used to obtain the loading on a cylinder 
following a cone. The reference diameter for cylinders is the cylinder diameter 
or D 0 = Dp = 2.17 meters. The cylinder of this configuration has a fineness 

ratio of 4.87. The curves of figure 5 are for cylinders following sharp cones; 
however, since the effect of blunting on the cylinder loading is secondary, 
these curves are also used for cylinders following blunt cones. From fig- 

for a cylinder following a cone with 

semi vertex angle Q n = 22.5° are obtained. Linear interpolation on 9 n is 

used to determine the appropriate values. These values are then multiplied by 
the cylinder diameter and tabulated in table B3* The final curve is shown in 
the following sketch: 


ure 5 (d.), the values of 


D ^ dx J 
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Cylinder Following Cone 

The frustum data are presented in figures 6 and 7* The parameters which 
affect the loading on a frustum are the nose angle 9 n , the preceding cylinder 
length the frustum angle (9f), and the Mach number (M). For the con- 

figuration under consideration, 9 n = 22.5, f a,l = ^* 87 , and 9-f = 5. The 
reference diameter D 0 for a frustum is the minimum diameter; therefore, 

D 0 = Dj_ = 2.17 meters. The fineness ratio of a frustum f^ and the local diam- 
eter D in terms of ' x/D stations can be obtained by using the formulas found 
in table B4 as follows, hence: 


f f = 2.42 


D = 


2.17(1 + 0.1750 



By using the foregoing equation the appropriate values of D can be calculated. 
Since f a is greater than 4.0, the loading on the frustum is no longer con- 
sidered to be a function of 9 n . Therefore, the only parameter for which 
interpolation of the available data is required is Mach number. The data of 
figure 6(g) and figure 6(i) were used with f a ^ = 4.0 and 9f = 5°- From 


these curves, the values of 



at the appropriate x/d q locations 


are read. These values are then multiplied by their respective local diameters 
to yield the final frustum loading curve. The desired data for Mach 2.4 are 
tabulated in the next to the last column on table B4. 
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Cylinder Following Frustum 


The data figures 8 (d) and 8 (e) are used to determine the loading on a cyl- 
inder following a frustum. For the configuration under consideration, the 
parameters affecting the cylinder loading are the frustum angle ( 9 f)> the 

cylinder-diameter ratio ^D^D Q and the Mach number (M) . Since curves are 

available for a 5 ° frustum angle, the only interpolations which are necessary 
are with respect to D^/D q and Mach number. The reference diameter for this 

component is 3*09 meters. The fineness ratio is 6 . 98 . The interpolated values 



are given in table B5- 


The normal-force-coefficient slope for the configuration shown in fig- 
ure B1 has now been determined for Mach 2.40. All that remains is to combine 
the component data. It should be noted that the values of x/D q at which the 
data curves were read are completely arbitrary. The only criterion to be con- 
sidered is the adequate description of the loading functions. When the com- 
ponent data are combined, discontinuities may exist in the final curve. This 
is a valid condition and these discontinuities are associated with the discon- 
tinuities which occur in the slope of a meridian section of the body. Where 
body-slope discontinuities do not exist, the loading function should be faired 
together to form a smooth curve. For this configuration this condition occurs 
at the sphere-cone interface. The reference point for plotting the final curve 
is also arbitrary, and, for this case, the nose of the vehicle has been chosen 
as the zero station for the p coordinate of the combined system. To convert 
the component x/D 0 coordinates to the combined coordinates p, the component 
x/d 0 values are multiplied by the appropriate D Q for each component and 
added to the final p station of the previous component. The initial and 
final stations for each component are given in the p columns of tables B1 
to B 5 . As stated previously, the reference area term is also arbitrary and 


should it be desired to convert the final 


dC 


N, 


f s ^£ 

l 


curve given in figure B1 


to 


a 


dx 


form, all that is required is division by the selected reference area. 


This may be the maximum cross-sectional area of the cylinder or the configura- 
tion planform area, or in the event fins are used, the fin planform area. 
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TABLE Bl.- APPLICATION OF METHOD TO SPHERICAL SEGMENTS 
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TABLE B 2 . - APPLICATION OF METHOD TO CONES 



Ho = 0.268 meter 


meters 

x - M-o 

Do D 0 

D, 

meters 

C dC N 

S< — — / , meters /radian 
at M = 2 . 14 - 

0.268 

0 

0.802 

I.67 

1.920 

.761 

2.17 

1+.51 
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TABLE B4.- APPLICATION OF METHOD TO FRUSTUMS FOLLOWING CONE -CYLINDERS 


D 0 = D^ = 2.17 meters 

e f = 5° 
e n - 22.5° 
f a,l * ‘‘•875 
d 2 ~ 3*°9 meters 


f = d 2 ~ P 1 = _L = 

f 2Dj tan 0^. D^ 

D = D 0 /l + 2 — tan 6 f ' 
V D n 


(i 0 = 12.499 meters 


2.423 






*( s 

ac NA 

dx f * 

{ 

’ dC No\ 

-j’ 




x M- “ M-o 

16 

J 




meters 

P° 

ll 

l«° 

meters 

l/radian at - 

meters / radian 

at 

- 




M = 2.18 

M = 2.81 

M = 2.18 

M = 2.4 

M 

= 2.81 

12.499 

0 

2.17 

0.23 

0.20 

0.50 

0.48 


0.43 

13.584 

.5 

2.36 

•57 

.12 

.87 

.66 


.28 

14.669 

1.0 

2.55 

• 50 

•15 

1.28 

• 95 


.55 

15.754 

1.5 

2.74 

.60 

•17 

1.64 

1.23 


.47 

16.839 

2.0 

2.95 

.66 

• 23 

1.93 

1.49 


.67 

17-757 

2.423 

5.09 

• 70 

•51 

2.16 

1.74 


.96 
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Figure Bl.- Product of reference area and the distributed normal-force -coefficient slope for 

typical launch -vehicle configuration. 
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table I 


COMPONENTS FOR WHICH LOADING FUNCTIONS ARE AVAILABLE AND 


THE PARAMETER RANGES COVERED 







CYLINDER FOLLOWING 
CONE OR OGIVE 


FRUSTUM FOLLOWING 
CONE- CYLINDER 



D 

o 





f* — , 1 — * — x 


BOATTAILS FOLLOWING LONG AND SHORT CYLINDERS 





Figure 1.- Parameters used to define geometrical components. 
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(a) M = 0.25 to 1. 

Figure 2.- Loading functions for spherical segments for various Mach numbers. 
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(k) Cylinders following 


Figure 5.- 



; M = 3-0 and 4.24 
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(e) m = 1.5; e f = io°. 

Figure 9 *- Continued. 
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